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One-dimensional chirality: strong optical activity in epsilon-near-zero metamaterials 
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We suggest that electromagnetic chirality, generally displayed by 3D or 2D complex chiral struc¬ 
tures, can occur in ID patterned composites whose components are achiral. This feature is highly 
unexpected in a ID system which is geometrically achiral since its mirror image can always be super¬ 
posed onto it by a 180 deg rotation. We analytically evaluate from first principles the bi-anisotropic 
response of multilayered metamaterials and we show that the chiral tensor is not vanishing if the 
system is geometrically one-dimensional chiral, i.e. its mirror image can not be superposed onto it 
by using translations without resorting to rotations. As a signature of ID chirality, we show that 
ID chiral metamaterials support optical activity and we prove that this phenomenon undergoes 
a dramatic non-resonant enhancement in the epsilon-near-zero regime where the magneto-electric 
coupling can become dominant in the constitutive relations. 

PACS numbers: 78.67.Pt, 81.05.Zx, 78.20.Ek 


The term chirality is generally used to express the 
asymmetric property of a 3D object which can not be 
superimposed onto its mirror image by translations or ro¬ 
tations. 3D chirality is an important feature in many or¬ 
ganic molecules (for example, 19 of the 20 common amino 
acids that form proteins are chiral) and its associated 
phenomena have enormous impact in several branches of 
science encompassing molecular biology, life science, op¬ 
tics, crystallography and particle physics. The concept of 
2D chirality also exists and a planar object is said chiral 
if it cannot be superposed onto its mirror image unless 
it is lifted from the plane. Although the 3D chirality is 
widespread in nature, examples of 2D chiral structures 
are very few jjlj. 


In the context of metamaterial science, artificial chi¬ 
ral structures, whose underlying constituents exhibit in¬ 
trinsic chiral asymmetry, have been investigated theo¬ 
retically and experimentally by several groups in 3D 
[g] and 2D configurations [6l4lQ|. Chiral metamaterials 
have attracted a good deal of attention since they can 
yield strong chiral bi-anisotropic response (due to cross¬ 
coupling between the magnetic polarization and the elec¬ 
tric one), giant optical activity, asymmetric transmission 
a repulsive Casimir force [12} and negative refractive 
index [ 131416 }. It is worth noting that electromagnetic 
chirality can also be present in the case where underly¬ 
ing constituents are not intrinsically chiral. In this case, 
the magneto-electric coupling results from the geometric 
chirality of the whole structure and the effect is driven 
by the radiation wave vector contributing to the overall 
chiral asymmetry (extrinsic chiralilty). 2D extrinsic chi¬ 
ral metamaterials have been first proposed by E. Plum 
et al. E2 and, in the considered media, authors have 
observed large optical activity which is indistinguishable 
from that occuring in media whose constituents are in¬ 
trinsically chiral. 


In this Letter, we show that ID systems can support 
a peculiar reciprocal bi-anisotropic response in the long 
wavelength regime. We evaluate the chiral tensor of a 
class of multilayered metamaterials whose constituents 
are intrinsically achiral and whose ID pattering is not 
characterized either by the standard 2D or 3D chiral 
asymmetry. We obtain the intriguing result that the 
chirality tensor does not vanish if the structure exhibits 
one-dimensional chirality , i.e. its mirror image can not 
be superposed onto it by using translations without re¬ 
sorting to rotations. Recently, a tunable bi-anisotropic 
response has been predicted in stratified media hosting 
graphene sheets [18}. In the general situation that we 
are considering here, the effect of the ID chiral asymme¬ 
try along the layers stacking direction is so pronounced 
that a medium with a three-layered unit cell, which is the 
simplest ID chiral structure, can be characterized by a 
marked first order spatial dispersion which mimics the bi- 
anisotropic response of a uniaxial omega medium jl9| . In 
addition, we analytically prove that ID chiral metamate¬ 
rials show optical activity (circular birefringence and cir¬ 
cular dichroism) and this phenomenon is enhanced in the 
epsilon-near-zero regime to the point to be as strong as in 
2D and 3D chiral metamaterials. Such non-resonant en¬ 
hancement is a consequence of the general fact that close 
to permittivity crossing points the nonlocal and possi¬ 
ble nonlinear polarization contributions can be compa¬ 
rable or even greater then the local linear part of the 
dielectric response. For example, a nonlinear epsilon- 
near-zero metamaterial can support efficient nonlinear 
processes [20|, Hl[ and novel classes of solitons [Hj. On 
the other hand, in linear epsilon-near-zero metamaterials, 
the medium response can be strongly affected by nonlo¬ 
cal response and it can be exploited to boost nonlocal 
phenomena such as the excitation of additional waves 
[23}. Here we suggest that nonlocal effects due to ID 
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FIG. 1: (Color on-line) Sketch of metamaterial slab (with 
N = 3 layers) and waves scattering geometry. The trans¬ 
mission amplitudes of right-handed (RCP,+) and left-handed 
(LCP,-) circular polarized plane waves are not equal for 6 0 

as manifestation of ID chirality. 


chirality are enhanced by the epsilon-near-zero condition 
without any resonant mechanisms (i.e. with no cavity or 
plasmonic effects). 

Let us consider electromagnetic waves propagating in 
a ID chiral metamaterial whose underlying structure has 
a unit cell obtained by stacking along the x-axis, N-layers 
of different media of thicknesses dj (j = 1,..., N) (see 
Fig.l). The electromagnetic field amplitudes E, H, asso¬ 
ciated with monochromatic waves (the time dependence 
e~ lut has been assumed where ui is the angular frequency) 
satisfy the Maxwell’s equations V x E = iui/i oH and 
V x H = —iwD, where the constitutive relations are D = 
e 0 eE, H = (l//r 0 )B and the relative dielectric permittiv- 

N 

ity e is a periodic function of period d = ]T) dj. In order 

j = 1 

to obtain an effective medium description of the electro¬ 
magnetic propagation valid in the regime where the ratio 
between the period d and the wavelength A is small, we 
exploit a standard multi-scale technique suitable for very 

Accordingly, 


general dielectric periodic profiles 2 
we introduce the parameter 77 = d/X -C 1 and the fast co¬ 
ordinate X = x /77 and, aimed at isolating the slowly and 
rapidly varying contributions, we consider the Fourier 
series of the dielectric permittivity and its inverse, i.e. 
e = (e) + E a n e ink ° x and e " 1 = (e" 1 ) + £ b n e ink ° x 

n^O njtO 

(where (/) is the mean value of the function / and 
ko = 2 'k/X). The basic Ansatz of our approach is 
A (r, X) = A(r) + SA(r,X) + r 1 [A^ (r) + <5A (1 ) (r,X)], 
where A and 8A are the slowly (averaged) and rapidly 
varying part of each electromagnetic field (E or H) 
whereas the superscript ( 1 ) labels the first order correc¬ 
tions in 77 . Substituting the Fourier series of e and e _1 
and the basic Ansatz into Maxwell equations, after sep¬ 
arating the slowly and rapidly varying contributions and 
retaining the terms containing 77 up to the first order, 
we obtain the equations describing the dynamics of the 
slowly varying parts of the electromagnetic field, namely 
V x E = iui/i qH and V x H = — 


D x — eo 


£|| Ex + ~j~(,9yEy + d z E z ) 
Kc 


Dy - CQ I Cj_Ey dyE x I , 


ft 

ko 


D z = e 0 ej _E Z - — d z E a 
ko 


(1) 


H = (l///o)B, e|| = (e *) 1 and ej_ = (e) and we have 
introduced the chiral parameter 


n^O 


1—nbri 


( 2 ) 


It is evident that in the limit 77 —>- 0 the chiral parameter 
ft vanishes and the multiscale approach considered in this 
paper reproduces the results of the well known standard 
effective medium theory. 

The obtained expression for ft of Eq. (121) highlights the 
key role played by the fundamental concept of ID chi¬ 
rality. A structure is ID chiral if its mirror image can 
not be superposed onto it by using translations with¬ 
out resorting rotations. From Eq.© we have that the 
parameter ft vanishes if the structure is ID achiral. In 
fact, in this case, the permittivity profile after a reflec¬ 
tion (say through the plane X = 0) can be superposed 
to the original dielectric profile, i.e. it exists a trans¬ 
lation X' = X + h such that e(X) = e(—X + h). It is 
straightforward proving that the corresponding dielectric 
Fourier coefficients are such that a_ n = e~xjp(—inkoh)a n 
and 6 _ n = exp(—inkoh)b n so that a- n b n = a n b- n and 
the series of Eq.© provides a vanishing ft. Therefore, the 
slowly varying and leading electromagnetic field can ex¬ 
perience the effect of the novel terms proportional to its 
first order derivatives in the effective constitutive equa¬ 
tions of Eq. © only if the multilayer does exhibit ID 
chiral asymmetry. 

Note that Eqs. m are of the kind Di = 
e 0 (tijEj + (XijndnEj) (i = x,y,z) describing media 
with a weakly spatial nonlocal dielectric response 
stemming from spatial dispersion [26|. In the present 
case the components of the tensor aij n are all vanishing 

except Oi Z xz — &xzz — ^/ ko , 0/.yxy — (Xxyy — ^/ ko 

and this agrees with the fact that a satisfies the anti¬ 
symmetric relation ttp = — aji n due to the Onsager 
symmetry principle. 

Exploiting the fact that the effective Maxwell’s equa¬ 
tions are invariant with respect to transformation D' = 
D - V x Q, and H' = H + iu> Q (where Q is an arbitrary 
vector) [23,|2a|, after setting Q = —eon/ko{E z e y — E y e z ), 
we obtain the equivalent effective constitutive relations 

D = e (e//) E--ft (e//)T H, 
c 

B = -ft (e//) E + /x (e//) H, (3) 


where = diag (e||, e± + ft 2 , e± + ft 2 ), = /i 0 / 

and 


we define 

(0 

0 

0 ^ 



ft( e// ) = 0 

0 


( 4 ) 
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or, in other words, the considered metamaterial is a re¬ 
ciprocal bi-anisotropic metamaterial and, specifically, it 
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FIG. 2: Frequency dispersion of effective electromagnetic pa¬ 
rameters. 


mimics the response of an omega medium since the chiral¬ 
ity tensor (reciprocal mag neto-electric coupling) is 

purely antisymmetric [19t |291 |. Standard chiral materials 
have, by definition, chirality tensor with not-vanishing 
trace and therefore we again stress that, in this Letter, 
one-dimensional chirality is purely geometrical feature 
of structures whose mirror image can not be superposed 
onto them by using translations without resorting to ro¬ 
tations. It is worth noting that the structure of the chiral 
tensor of Eq. (J3J) of a one-dimensional chiral metamaterial 
can be directly obtained from symmetry considerations 


The stratified medium whose unit cell has two layers 
(TV = 2) has vanishing chiral tensor, n = 0, since it 
is always possible to superpose the mirror image of 
the structure with the original one through a suitable 
translation. Accordingly, as discussed in Ref. [33}, 
standard metal-dielectric stratified metamaterials show 
purely second order spatial dispersion. Therefore, the 
simplest stratified medium which can show not vanishing 
chiral tensor has three different layers in its unit cell 
with different permittivities 61 , 62 , 63 . The Fourier 
coefficient of the dielectric permittivity e(X) are a n = 
2 k H £ i - £ 3 ) + (d - e 2 )e-«i + (e 2 - e 3 )e‘ 2 ™*], 
where n ^ 0 and /, = di/d is the filling fraction of 
the i-tli layer, and those of l/e(X), b n , have the same 
expression with the three permittivities 6 j replaced by 
their inverses 1/e^. Inserting these coefficients into 
Eq.® and summing the resultant series using the re¬ 
lation J27=i sin{2 ™ fi) = x ( 2 f? _ 3 /* 2 + fi) we obtain 

«=tt^ii (ei - £2)( : 1 i ~: 3 3 )(£2 - e3) /i/2/ 3 . 

As an example, we consider the electromagnetic ef¬ 
fective response of a periodic three-layer structure with 
fi = 0-18, h = 0.33, h = (1 - fi ~ /a), d = 25 nm, 
where we choose the slab layers of type 1, 2, 3 to be filled 
with Ag (1), Si0 2 (2), PMMA (3), respectively. The sil¬ 
ver dielectric permittivity is given by the Drude model 
6 i = 6 b — w 2 /(w 2 + iauj) (where ej, = 5.26, a = 7.06 • 10 13 
s _1 , uj p = 1.45 • 10 16 s' - 1 0 ), whereas the frequency 
dispersion of Si0 2 (e 2 ) and PMMA (e 3 ) are described by 
empirical Sellmeier equations with different coefficients, 
respectively EHH- In Fig.2, we report the effective di¬ 
electric permittivities (panel (a) and (b)) and the chiral 
parameter k (panel (c)) as function of the wavelength A. 


As shown in Fig. 2(a), the ID material displays an ef¬ 
fective dielectric permittivity with a zero-crossing point 
near A ~ 0.51 p.m (Re(ej_) ~ 0) and it is characterized 
by a hyperbolic response (Re(e±) < 0, Re(t\\) > 0) for 
A > 0.51 pan. 

The discussed electromagnetic bi-anisotropic response 
pertaining ID chiral metamaterials is able to support 
specific optical activity effects such as circular birefrin¬ 
gence and dichroism. Besides, the permittivity ex of the 
considered stratified structure, as above discussed, can 
have zero-crossing points (if both metal and dielectric 
layers are used) and therefore the epsilon-near-zero 
regime, i.e. |ej_| -C 1, can occur in a spectral region 
around the crossing-points [ 3 TI 132} where non-resonant 
enhancement of the chiral response is expected. In 
order to appreciate the effects of the ID chirality 
and their enhancement due to the epsilon-near-zero 
regime, we consider a slab of ID chiral metamaterial 
of thickness L illuminated by circularly polarized plane 
waves impinging from vacuum with incidence angle 6 
(see Fig.l). We here do not consider a substrate which, 
although important for practical implementations, does 
not provide an essential physical ingredient in the 
mechanism discussed here. For a given transversal wave 
vector k y , there are two forward eigenwaves excited 
within the slab, i.e. the ordinary and extraordinary 
waves E(°) = (k { z o) e y - k y e z ) e ik ^ z + ik vV , E( e ) = 


Eie) 

kyGy “f- kz 

jV. (e) 

z+iky-y 

5 

where 

k = ijk%e±- 

fc 2 and hi'^ ^ 

l k o e \\ 

( 1 + S) 



By imposing the continuity of tangential electric E 
and magnetic H 7 fields at the vacuum-metamaterial 
interfaces, one can obtain the complex transmission 
matrix t relating the transmitted E^) and incident 
E^) circularly polarized fields, namely = ti m Em , 
where indices m,l correspond to polarization states 
of transmitted and incident wave, which can be ei¬ 
ther right (+) or left (-) circular polarizations. To 
quantify the optical activity of the considered meta¬ 
material slab, we evaluate the difference in magnitude 

A = 2(|i_| 2 — |t ++ | 2 )/(|t_| 2 + |t++| 2 ) and in phase 

<5$ = arg[t _) — arg(t ++ ) of the two diagonal terms 

of the transmission matrix, since these two paramaters 
are associated to circular dichroism and circular bire¬ 
fringence, respectively. Using the effective permittivities 
and chiral parameter reported in Fig.2 for a slab of 
thickness L = 0.4 /rm, in Fig.3 we plot the percentage 
difference A (panel (a)) and (panel(b)) as functions 
of the wavelength A for 9 = 0 deg and 9 = ±30 deg. 
It is evident from Fig.3 that close to A ~ 0.51 /im, 
which corresponds to a crossing point of the permittivity 
e± (see panel (a) of Fig.2), both circular dichroism 
and circular birefringence show marked peaks. From 
the first of Eqs. Q: it is worth noting that in the 
epsilon-near-zero regime |ej_| <C 1 the magneto-electric 
coupling terms containing the chiral parameter can be 
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FIG. 3: Optical activity of a slab whose permittivities and 
chiral parameter are reported in Fig.2 and of thickness L = 
0.4 pm. (a) Circular dichroism (A(%)) and (b) circular bire¬ 
fringence as functions of wavelength A for 9 = 0 deg (solid 
line), 9 = 30 deg (dash-dot line) and 9 = —30 deg (dash line), 
(c) Parameter p measuring the magnitude of the magneto¬ 
electric coupling term relative to the local dielectric contribu¬ 
tion from the //-component of the first of Eq.Q. (d) A(%) 
(squares) and (circles) as function of 6 for A = 0.5 /im. 



FIG. 4: Comparison between effective medium theory (EMT) 
(solid lines) and full-waves analysis (dashed lines) predictions 
about the optical activity of a metamaterial structure whose 
material parameters are those used for obtaining Fig.2 with 
L = 0.4 pm and 9 = 30 deg. 


comparable or greater then the local dielectric contribu¬ 
tions and this triggers the non-resonant enhancement of 
the effects resulting from ID chirality. This mechanism 
is proved in Fig.3(c) where we plot the parameter 
p = \eoK 2 E y + i^H' z | / \eoej_Ey |, measuring the magni¬ 
tude of the magneto-electric coupling term relative to 
the local dielectric contribution from the //-component 
of the first of Eq.(|3]) and evaluated at z = L~ in the 
same situation of Fig.2(a-b), and p shows a marked peak 
around A = 0.51 /mi reaching the maximum of 1.2. The 
same mechanism can also be grasped from the expression 
of the extraordinary eigenwaves where it is evident that 
the first order parameter ruling the effect of ID chirality 


on the electromagnetic field is n/e±_ rather than k and 
this produces the enhancement of optical activity in 
the epsilon-near-zero regime. In Fig.3(d), we plot the 
percentage difference A and (5$ as functions of the angle 
9 for A = 0.5 pm. Note that the circular dichroism 
and birefringence are completely absent for 9 = 0 deg 
(black solid line in Fig.3(a-b)) and they present opposite 
signs for opposite angles (dash-dot and dash lines in 
Fig.3(a-b)). These features are due to the fact that 
the //-component of the extraordinary wave (absent 
in local uniaxial materials) is proportional to k y . ID 
chiral metamaterials share the same angular dependence 
observed in 2D extrinsic chiral metamaterials E2 and 
they show optical activity analogous to one occurring in 
2D and 3D chiral materials. 


In order to check the validity of the above results 
we have performed 3D full-wave simulations 37] where 
two circularly polarized monochromatic plane waves are 
made to interact with the metamaterial structure as 
shown in Fig.l and optical activity is accordingly re¬ 
trieved. The considered structure has the same material 
parameters as those used to obtain Fig.2 and the scatter¬ 
ing process has been numerically investigated for L = 0.4 
pm and 9 = 30 deg. In Fig.4 the comparison between the 
results of the effective medium theory (EMT) of Eqs.([3]) 
and those of the numerical simulations are plotted and 
their good agreement is evident, together with the ENZ 
enhancement of the optical activity. The only slight dis¬ 
crepancies appear in a narrow spectral range centered at 
the ENZ wavelength A = 0.51 pm since, in such spectral 
range, the effects of the (here neglected) higher order con¬ 
tributions to the dielectric response (e.g. containing sec¬ 
ond order electric field derivatives in Eqs. dT|) ) are boosted 
(as shown in Ref. II) as well as those due to first order 
contributions discussed in the present letter. In addition, 
such discrepancies can also arise from using an homoge¬ 
nization theory performed for an unbounded medium in 
the presence of a slab which has boundaries. 


In conclusion, we have shown that a multilayered struc¬ 
ture exhibiting ID chiral asymmetry along the stacking 
direction shows a reciprocal bi-anisotropic electromag¬ 
netic response whose effects are comparable to those as¬ 
sociated with 2D and 3D complex omega media. Here 
the bi-anisotropic response is supported by ID metama¬ 
terials which are much easier fabricated and theoretically 
investigated than 2D and 3D chiral metamaterials. In ad¬ 
dition, we have proposed a strategy to enhance the op¬ 
tical activity stemming from ID chirality in metamateri¬ 
als by exploiting the epsilon-near-zero regime where the 
magneto-electric coupling contribution to the medium di¬ 
electric displacement is even greater then the standard 
dielectric permittivity one. Our findings could be essen¬ 
tial for conceiving ultra-efficient bi-anisotropic photonic 
devices in the optical frequency range. 
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